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_|- , We consider the creation of non-zero Chern-Simons number in a model of the early 

• Universe, where the Higgs field experiences a fast quench at the end of inflation. 

We perform numerical lattice simulations in the Abelian Higgs model in 1+1 di- 
mensions and in the SU(2)-Higgs model in 3+1 dimensions with an added effective 
, CP-violating term. We also comment on the appropriate choice of vacuum initial 

O ; 

I ' We consider a scenario of baryogenesis in the early Universe, where 

inflation ends at the electroweak scale 1 ' 2 > 3 . After inflation, the fields are 
, in the ground state at zero temperature (no reheating from the inflaton). 

Higgs symmetry breaking happens through a quench causing the Higgs 
field to go through a spinodal instability. CP (and C) is broken by an 
effective term in the action. Baryon number (B) is not conserved due to 

■ Chern-Simons number (Nqs) changing transitions in the electroweak sector, 
through the relation B(t) - B(0) = 3 [N cs (t) - iVcs(O)]. The three criteria 
for baryogenesis are thus fulfilled. 

We study two models of this system: The Abelian-Higgs model with 
C and P violation in 1+1 dimensions and the SU(2)-Higgs model with 
CP violation in 3+1 dimensions. We apply the classical approximation 
in which quantum averages are replaced by averaging over an ensemble of 
classical initial conditions, evolved using classical equations of motion. This 
approximation is valid since large occupation numbers Uk are generated by 
the spinodal instability (see below). 

We assume that the end of inflation triggers electroweak symmetry 
breaking by a rapid change of the effective mass in the Higgs potential 
V{<j>) — + A(cfy (f) 2 . We model it by an instantaneous change of 
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the sign of the mass at t = 0, — > — /i 2 . For i > the original ground 
state is unstable, and momentum modes with |fc| < [i grow exponentially 
through the spinodal instability. This proceeds until the quartic term be- 
comes important and the fields eventually thermalize in the broken phase 
at a low temperature. 

For small couplings, we can initially neglect the quartic term and gauge 
fields and solve the quantum evolution of the unstable Higgs field exactly 
4,5 . In terms of correlation functions, the result is (ujf — ^/k 2 ± fx 2 ): 




where we have defined effective time-dependent particles numbers nt, hk 
and frequencies Wfc. For k in the unstable region (\k\ < fi), = i-yV 2 — ^ 2 
and nk and hk grow exponentially cx exp ^2-\//i 2 — k 2 tj . At some time 
during the rolling down, we can switch to a classical description through an 
ensemble of classical initial conditions. The classical correlation functions 
are taken to be the real parts of the quantum ones. 

The initial quantum state (t = 0) has nk = hk = 0. For t > 0, hk ^ 
and the uncorrelated variables are £^ = (ir-k/^/uJk ± y^LUkfik) /V2. At t = 
t ro we draw sets of random numbers £^ from the distribution: 

As t grows, fik — > rik + |, and we obtain the classical relation ^ — > 0, or 

The classical approximation of the dynamics is valid for sufficiently large 
t ro - However, if the assumption of small coupling is correct, the classical 
evolution is equal to the quantum spinodal evolution, and it should not 
matter which t ro we choose. A special case is t ro — 0, where nk + \ = \- 
We call this the "Just the half" case. For all t ro , we only initialize modes 
with |fc| < /i. Eventually, the quartic term will become important and the 
approximation breaks down. t ro can therefore not be too large (Fig. 1, left 
plot). In the bulk of our study of the CP- violation we chose to describe 
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Figure 1. Left: (<j>*<j>) for different values of t ro . At the larger roll-off times, the behavior 
is significantly different from the "Just the half" method, indicating that the quadratic 
approximation breaks down. The coupling is A/^t 2 = 1/8. Right: (Ncs) for different 
initial condition schemes. The thermal scheme is quantitatively different, but still shows 
the important features, including the "initial dip" . The spinodal case at later and later 
t ro differs in the initial dip due to the coupling to the C and P violation being introduced 
later (at t ro )- 



the state just before the quench by a Bose-Einstein (BE) distribution at 
low temperature: 
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The Abelian Higgs model in 1+1 dimensions is defined by the action: 



S = 



dtdx 



1 



Y^F^F^ + {D^YD^ - M 2 0*</> + \(<t>*4,y 



4e 



The last term violates C and P. We used BE initial conditions at T /ran = 
0.1 with A/^t 2 = 1/8. The initial gauge field is given by solving the Gauss 
constraint. We used averages over an ensemble of 1000 initial configura- 
tions. 

For non-zero k we can estimate the initial effect of the C/P break- 
ing term from the equations of motion, with the result Ncs(Un) = 
K(m 2 v Lm_y)/(m| / 6v / 87r), where {4>*4>) (U n ) = [i 2 /3\. Indeed we see such 
an initial dip (Fig. 1, right plot). The final (at tmn = 600) (iVcs) is linear 
in k (Fig. 2) and depends sensitively on mn/mw (Figs. 2,3). 

The SU(2)-Higgs model is given by the action: 



Sri — — 



d A x 



+ A 



-K^tj/TrF^F, 
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Figure 2. Left: The dependence of < Ncs >final on K ls unear (here for mu/ray/ = 
1. Right: The dependence on the mass ratio of the Higgs and W particles is very 
complicated. 




Figure 3. Depending on the Higgs to W mass ratio, oscillation frequencies seem to 
conspire to give smaller or larger final iVcs; even to the point of giving opposite signs. 
Here mu/mw = 0.63 (left) and 1.06 (right); n = —0.05. 



F^v is the dual field strength and 4> the Higgs doublet field. The last term 
breaks CP. We used initial conditions similar to the 1+1 dimensional case. 
The results for ran = rnw are shown in Figs. 4, 5. We can estimate the 

2 4 

final temperature after symmetry breaking from ^g^T 4 = V(0) = ygj. 
Using n B = s = ^-g^T 3 with g* = 80, Lm H = 21, A = 1/18, we 

find (k = 167t 2 k) 

< Ncs > k =3/ m 2 ^ 0.2 =*> — (k) ~ 7 x l(T 3 KmL ~ 5 x l(T 5 KTeV 2 

' H rij 

This result seems quite reasonable, and requires that k ~ 10 _7 /m^. For 
additional information on tachyonic preheating, see the accompanying con- 
tribution 6 . 
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k=3/m H , Ncs and Winding. 
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k=0, Ncs and Winding. 



Figure 4. Two examples of trajectories in the SU(2)-Higgs simulation with k = 16tt 2 k = 
and k = 3/mjj. Left: For this initial condition the effect of the CP- violation is a small 
shift in Ncs and winding number. Right: Here, the trajectories of both observables are 
split to end up at different integer values. 
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Figure 5. Left: Distribution of Higgs winding number and Chern-Simons number for 
k = Right: Bootstrapping the average Chern-Simons number (^Vcs) nna i- Notice 

that since we used the same initial conditions for zero and non-zero k, the effect of CP- 
violation is the difference between the two Gaussians. Due to statistics, these are both 
shifted to the left. 
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